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Generic properties of elastic phonon transport at a disordered interface 

are studied. The results show that phonon transmittance is a strong function 
of frequency and the disorder correlation length. At frequencies lower than 
the van Hove singularity the transmittance at a given frequency increases as 
the correlation length decreases. At low frequencies, this is reflected by dif- 
ferent power-laws for phonon conductance across correlated and uncor related 
disordered interfaces which are in approximate agreement with perturbation 
theory of an elastic continuum. These results can be understood in terms of 
simple mosaic and two-colour models of the interface. 

I. INTRODUCTION 

The transfer of phonons across the interface between two solids is crucial to the operation 
of a range of novel devices from semiconductor nanostructures, quantum-cascade lasers, 
vertical structures, superlattices, quantum dot arrays to the cryogenic phonon-mediated 
detectors of elementary particles. In such structures non-radiative transitions and hot- 
electron thermalization gives rise to fluxes of highly energetic phonons spanning the whole 
of the Brillouin zone. 

In early studies |IJ , phonons were regarded as bulk elastic continuum waves so that trans- 
mission and reflection of energy at the interface between two materials were a consequence 
of acoustic mismatch. Such a macroscopic continuum theory describes generic low-frequency 
properties but does not account for important features such as phonon dispersion, the Debye 
cutoff and weak bonding at the interface, for which a lattice-dynamical model is required 



0. Most early investigations assume that the solid-solid interface possesses in-plane transla- 
tional invariance, whereas, it is known that disorder-induced scattering can play a significant 
role 0>f§. Within a lowest order perturbation theory ||, it is predicted that disorder can 
induce a strong frequency dependence in phonon transport coefficients with a power-law 
crossover determined by interface conditions. The aim of this paper is to go beyond this 
elastic continuum perturbative approach, by formulating an exact lattice-description valid 
for terahertz phonon transport, which includes features such as phonon dispersion and sin- 
gularities in the phonon density of states. 

To formulate the problem of phonon transmission (reflection) we consider two identi- 
cal semi- infinite leads of uniform cross-section attached to a 'scattering region'. The two 
semi-infinite perfect crystals are treated as waveguides for incident and scattered phonons 
whose dynamics is described by the lattice equations of motion. The geometry analyzed 
is shown in Fig.l, where the scattering region consists of a single disordered atomic plane. 
Disorder at the planar interface is introduced microscopically as correlated or uncorrelated 
random variation of the on-site atomic masses and we assume that there are no significant 
anharmonic effects. To calculate transport properties of such a system, we use methods 
specifically developed to investigate electronic transport through phase-coherent structures 
f^H- Transport coefficients are calculated as a function of phonon frequency within the 
Landauer-Buttiker formalism |J by employing an exact recursive Green function technique 

fi- 
ll. DETAILS OF THE CALCULATION 

In what follows we present the results for two lattice dynamical models: an fee lattice 
with central-force nearest-neighbour interactions and a scalar square lattice. For the former 
the atomic masses M/ mn are located at equilibrium positions r; mn = (la,ma,na), where 

1, to, n sum to an even number and a\/2 is the nearest-neighbour spacing, whereas for the 
square lattice a is the lattice constant and n = 0. In the latter case we do not distinguish 
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between phonon polarizations and hence, are not concerned about mode conversion due to 
scattering. 

In the harmonic approximation the linearized equations of motion of an atom located at 

r /mn IS 

Ml mn itf mn = — ^ ^lmn,l'm'n' U Vm'n'^ (1) 
P, I'm'n' 

where uf mn is the atomic displacement vector, 
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and Yvm'n' — ?imn is taken from the first coordination sphere. The indices a, (3 label the 
x, y, z components of displacements for the fee model and should be dropped for the scalar 
model. 

All atoms in the semi-infinite leads are of unit mass, whereas the random masses at 
the interface have mean value unity and standard deviation a. For simplicity we allow the 
masses to fluctuate only along a single direction, which is chosen to be the [101] and [100] 
cubic crystallographic directions for the fee and square lattice respectively. Henceforth, we 
refer to this as the 'disorder- line' and denote masses at different positions along the line by 
Mj. Thus, even for a three dimensional lattice, due to the residual translational invariance, 
scattering is confined in the plane fixed by the disorder-line and the angle with respect to 
the interface determined by the conserved component of the incident phonon momentum. 
From this point of view our scattering problem is essentialy two dimensional. 

An ensemble of iV positive masses, where N is the number of atoms along the direction 
of the disorder-line in the lead cross-section, is generated as follows: First we introduce the 
correlated random numbers 

Xj = J2( a k cos(kj) + b k sm(kj)), (3) 

k 

where a k and b k are Gaussian random numbers with zero mean value and (ak(ik>) = 
(bkbk') = $kk'Z ~- k is taken commensurate with the phonon wavenumbers determined 



by the boundary conditions in the transverse direction and £ denotes the correlation length. 
Random masses Mj along the disorder-line with mean value unity and standard deviation 
a are obtained from %j D Y the linear transformation 

Mj = 1 + a ■ -MZML . (4) 

In order to isolate the effect of correlations, we compare phonon transmission (reflection) 
coefficients for such correlated interfacial disorder, with those of an uncorrelated disorder 
configuration obtained by randomly 'shuffling' the above set. 



III. S-MATRIX AND PHONON INTERFACIAL TRANSMITTANCES 

In the absence of inelastic scattering, phonon transport through an arbitrary scattering 
region can be described by the scattering matrix S ||, which depends on the conserved 
phonon frequency and yields probabilities for all scattering transitions within a system de- 
scribed by a dynamical matrix D^ n l , m , n , = M^K^ Vm , n ,MySn<- If the scattering region 
is connected to external reservoirs by phonon waveguides with open channels labelled by 
quantum numbers u, then S'-matrix elements S uu /(uj,D) are defined such that \S ul/ /(u, D)\ 2 
is the outgoing flux of phonons along channel u, arising from a unit flux incident along 
channel z/. The scattering matrix defined this way is unitary, i.e., SS' = 1, and since 
S uu >(u,D) = S u i u (lu, Z?t), where D is real, S(u,D) is symmetric. If u, i/ belong to the left 
(right) lead then the corresponding S'-matrix elements are left (right) reflection amplitudes 
which we denote r vv i = S vv i{r' vv , = S uu >), whereas, if u, i/ belong to the right (left) and 
left (right) lead, the S'-matrix elements are left (right) transmission amplitudes denoted 
t V v' — S vv i(t' w , = S vu i). It is evident that for an iV'-mode waveguide three distinct trans- 
mittances can be defined, 

N' N' 

T vu l = \t vu '\ , Tyl = \tuu'\ ) T = \tvv'\ ) (5) 

V=\ v ,v' = \ 

To isolate the frequency dependence of the overall phonon transmittance it is natural 
to normalize T by the density of states. In the presence of finite-width leads this means 
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averaging T over all open scattering channels N' and is equivalent to the angle-averaged 
transmission coefficient, 



(T(lu,D)) = T(lu,D)/N\ 



which is related to the corresponding reflection coefficient by 



(6) 



(R{u,D)) = l-(T{u,D)). 



(7) 



In what follows the S'-matrix is obtained by solving Dyson's equation for the Green 
function at the surfaces of the two leads (see Fig.l for definition of L,R), namely 



' G(L,L;uj 2 ) G(L,R;uo 2 ) » 



G(R,L;u 2 ) G(R,R;u< 



,2\\-l 



G{u 2 ) = {g{u 2 )^-D eff {u 2 )) 



where g{uj 2 ) is the surface Green function of the semi-infinite leads 



(8) 



(9) 



9[u 



g(L,L;u 2 ) 
g(R,R;u; 2 ) 



(10) 



and D e ff(uj 2 ) an effective dynamical matrix 



Deff(u; 2 ) 



D eff (L,L;uj 2 ) D eff (L,R; 



OJ 



(11) 



/ 



D eff (R,L;ou 2 ) D eff (R,R;u 2 ) 

The latter describes the effective (renormalized) coupling between the surfaces of the two 
leads, obtained by projecting out the internal degrees of freedom of the scatterer using the 
recursive Green function technique. The Green function g(u) 2 ), the spectrum and structure 
of the allowed wavequide modes at fixed phonon density of states, as well as S(u, D) are 
calculated using a general algorithm developed in Ref. M. 
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IV. PHONON TRANSMITTANCE AT LOW FREQUENCIES 



We present results for the average transmission coefficient (T(u)) c where the subscript 
c denotes a configurational averaging. In Fig. 2 (a), (b) the frequency dependence of the 
overall transmission coefficient (T(u)) c averaged over 10 realizations of the interfacial mass- 
disorder for the fee lattice and over 500 for the scalar square lattice is shown. The width of 
the disorder-line is 29 and 100 nearest-neighbour spacings respectively, which corresponds 
to 3 x 29 and 1 x 100 degrees of freedom. 

For the fee model periodic boundary conditions in the lateral direction (disorder-line) 
were used, whereas for the scalar model we adopted fixed-end boundary conditions. The 
boundary conditions set the precise frequencies at which propagating waveguide modes be- 
come available in the external leads, but apart from features which scale inversely with the 
number of such modes do not affect the shape of the transmittance curves. Figure 2 shows 
that for both fee and scalar lattices, phonon transmittances (reflectances) follow the same 
universal behaviour. At sufficiently low frequencies, the phonon transmittance for correlated 
disorder is lower than for uncorrelated interfacial disorder with the same mean mass and 
standard deviation a . 

At low frequencies the results of lattice-dynamical calculations in Fig.2(a),(b) are in 
good agreement with perturbation theory ||. To model the disorder in a thin layer of 
dissimilar material sandwiched between irregular boundaries was introduced. The density 
and elastic constants of this layer were assumed to differ from those of the bulk on either 
side of the layer. To make a comparison with the perturbative analysis we consider the case 
of no fluctuation in the elastic constants and a line of random masses at the location of 
the idealized interfacial plane z = 0, i.e., Ap(x, z) = Ap((x)6(z) = Ap(x)5(z), with £(x) a 
random function of lateral coordinate with zero mean value. Apart from lattice discreteness 
which plays no significant role in the long-wavelength limit there is no difference between 
our model of site-to-site mass disorder and the model of interfacial roughness of Ref. 0. 
The following relation is useful when comparing the results 
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M 2 (O 



(12) 



p ) « 2 ' 



For such a two dimensional model the perturbation theory || with uncorrelated disorder 
£ = a and 1, where k is the phonon wavevector, yields 



while for correlated disorder at sufficiently high frequencies so that k£ 3> 1 the pertur- 
bative result is 



More generally for all u, Fig. 2(c) shows the perturbative results for £ — ► and £ — ► oo. 
Equations (|T3"D and (|14]) are in good agreement with the exact lattice-dynamical calculation. 
They both show lower transmittance for correlated disorder over a wide range of frequencies 
until the faster power dependence of phonon transmittance for uncorrelated disorder causes 
the two curves to cross at a frequency uj* . According to the perturbative calculations u>* 
should be close to half the maximum phonon frequency u> max independent of disorder char- 
acteristics, correlation length and strength. Despite the similarities, this is not evident in 
the lattice-dynamical calculations. Indeed, the crossing point is quite close to a singularity 
in the phonon density of states and thus using an elastic continuum model at frequencies of 
order u* is an oversimplification. 

Another simplifying assumption of the theory in Ref. j5] is the neglect of multiple scat- 
tering. This gives the power law (R(u))^ nc ~ u d+l (Rayleigh scattering) with subsequent 
decrease by one (two) power (s) for correlated disorder in 2D (3D). A log- log plot of phonon 
reflectances showed that both lattice dispersion and multiple phonon scattering effects re- 
sult in a change of the power exponents to a slightly different value. However, the drop by 
approximately one power for our 2D scattering problem in the case of correlated disorder 
prevails in all cases. The effect is attributed to the restrictions in the phase volume available 
for the scattered states, due to the correlation-induced finite width of the disordered spectral 
distribution. 




(13) 




(14) 
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V. SPECTRAL PROPERTIES OF PHONON TRANSMITTANCE 



Phonon transmittances over the whole frequency spectrum for the scalar model for two 
different values of disorder strength a are shown in Fig. 3. In addition to the low frequency 
behaviour discussed in Sec. IV, the transmittance shows a cusp at a frequency corresponding 
to a van Hove singularity in the phonon density of states and rapidly drops to zero at the 
maximum frequency of the phonon spectrum u max . At the van Hove singularity the difference 
between correlated and uncorrelated disorder realizations is much less pronounced. All these 
features are also seen for the fee model but the exact behaviour at van Hove singularities 
is less clear due to the finite-size effects described in the previous section. To interpret the 
results we concentrate on the scalar model. 

The general expression for the transmission amplitude from incident channel v' to channel 
v for systems with nearest-neighbour interactions has been derived in ||. The explicit 
expression for the scalar model with periodic boundary conditions is 

1 N frT 

*»" = T-A E e-^G LR (m,m 1 )W(m h m , )e^ m '%^e^ a . (15) 
-<V , i V V,,i 

m,mi,m'=l V " 

Here v v is the component of the group velocity in x-direction along the lead axis of the 
z/-th channel mode, m labels sites at the interface (along lateral direction), and G^ R has 
been defined in Eq.(H). The matrix W is given by H 



2 N 

W{m 1: rri) = -= J2 e ikV ^ mi - m ' )a (e ik > - e~ ik >). (16) 



In formulae (T5]TE|), k y v = as a result of transverse quantization with periodic bound- 
ary conditions, while denotes the wavenumber along the lead axis for phonons of fixed 
frequency uo of state v (z/-th channel). 

Using the complete orthonormalized set of functions {-^e lk ^ a } describing the profile of 
modes in the transverse direction, GLR(ui,mi) can be written 



1 



N 



G LR (m, mi ) = ^ £ e* k > a G LR ^^')e-^' mia . (17) 



S 



By substitution of Eg. (]T6| ) ; ( ]T7|) to ([T5| ) one obtains 

2 I 

t w > = ^G LR (u, v') % S in{kl,a)J^e- ik >. (18) 

Thus, to calculate the channel-to-channel transmittance (T uu i) c = (\t uu >\ 2 ) c or generally 
the overall transmittance (T) c one needs to consider the non-trivial statistically averaged 
product of two Green functions, (G\l_r(V, v'}G* LR {y, z/)) c or (GL R (m,m')G* LR (s, s')) c . 

In what follows, we shall demonstrate that generic features can be captured by consid- 
ering the case of a regular array of atoms with mass M' located at the interface, thereby 
forming a linear (planar) defect in a two(three) dimensional system. 

In this case, the full Green function Glr{v, v') is of the form 

O M fr-0= * ( ??f' (19) 
1 - a'uj 2 g {is, v') 

where g Q (i/, z/) = —4i5 u y j 1 {oj^ nax sin{kli)) is the Green function of the infinite ideal lead 
at coinciding arguments, and a' = M '^ I M is the relative mass difference at the location of 
the defects. Substituting this result into Eq.(|(J) and ([X8|) yields for the transmittance 

(T) = — V sm^Ka) 

v(open) \ v i \ u) max I 

Eq.(|20|) exhibits the following features: (i)(T) equals unity for an ideal interface with 
a' = 0. (ii)At low frequencies both uj, —>■ 0, sin(k^a) ~ u and therefore (T) — > 1 — au) 2 , 
where a is numerical coefficient, in full agreement with the results of perturbation theory 
H and formula ([HP- (iii)At the end of the phonon spectrum (T) — > since propagating 
channels are closing with sin(k^a) approaching zero at the Brillouin zone boundaries. This 
effect is analogous to the total reflection by a single mass defect in one dimension. (iv)At 
the van Hove singularity all channels are open and Eq . (p0|) reduces to 

1 * ^n 2 (— ) 
(T) = A^?ism 2 (^)+4a' 2 ' (21) 

which suggests that for a 1 < 1 the dependence on disorder strength is much weaker 
than at low frequencies. Also, sin(k^a) is small for a large number of contributing channels 
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thereby yielding the characteristic cusp. (v)We also note that Eq.(|20|) has been derived for 
the [100] orientation of the interface but since the phonon group velocity along the lead axis 
is proportional to sin(k*a) the transmittance is in general anisotropic. 

In Fig. 4(a), we plot the phonon transmittance given by formula ( p0| ) for selected values of 
a' along with the numerical results for correlated disorder realizations. Notice that in spite 
of the differences in the type of interfacial defect the curves look very similar. A natural 
way to explain this result is to view the interface for correlated disorder as a mosaic of 
clusters of average diameter ~ £. Suppose that within a particular cluster i all sites are 
occupied by defects with the same mass and characteristic relative mass difference a\. As 
long as the phonon wavelength is smaller than the correlation length, i.e., 3> 1, the overall 
transmittance can be considered as the average over all clusters with local transmittances 
given by Eq. (p0|) . In Fig.4(b) such an average is plotted for clusters with excess mass fully 
balanced by clusters with mass deficit. The relative mass difference is taken from a Gaussian 
distribution with width equal to the disorder strength a. Clearly, this procedure gives a 
better fit to the correlated disorder curve, which suggests that the phonon transmittance 
across a correlated disordered interface can be split to a statistically averaged mean field 
contribution induced by mosaic grains and a fluctuating component. The latter is controlled 
by the rapid residual variations in the relative mass difference at the cluster boundaries as 
well as the fluctuations in cluster sizes. 

The qualitative agreement between correlated disorder realizations and the mosaic pic- 
ture prompts us to ask whether there exists a simple model for the case of uncorrelated 
disorder. To explore this possibility, we examine a two-colour disorder model where linear 
interfacial defects are built up from atoms with relative mass difference a' which alternates 
sign every L/a sites. This model is motivated by the observation that uncorrelated disorder 
contains all spatial harmonics with equal probability. In what follows, we compare the exact 
transmittance for uncorelated disorder with that of the two-colour model, averaged over 
many values of L. In Fig. 5, we plot the mean arithmetic value of the phonon transmittances 
for L = 2,4, .., 10a taken after averaging over 500 two-colour disorder realizations given by 
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a Gaussian distribution with zero mean and y (tf" /2 ) = o. The remarkable agreement with 
the calculations for the uncorrelated disorder over a wide range of frequencies suggests that 
the mean field contribution for this case can be calculated by picking up few elements from 
the basis set of the two-colour disorder model. 

Finally, we compare the results of our calculation with the diffuse mismatch model DMM 
introduced in || and extensively used in studies of Kapitza conductance [[3|,j9|,|IO] . Application 
of the DMM at a boundary between two solids with identical acoustic properties gives a 50% 
phonon transmittance independent of frequency. It can be easily demonstrated that the 
Kapitza conductance ok for our model is larger than that predicted by the DMM. It is also 
less than the radiation limit which in this case would coincide with acoustic mismatch theory 
(total transmission). Moreover, ok is larger for uncorrelated disorder than for correlated. 
Finally, one notes that the results of [|2| which are valid for phonon transmission through 
the ideal interface between two different fee crystals a transmittance which is frequency- 
independent over the whole spectrum up to a cut-off frequency corresponding to maximum 
phonon frequency in one of the crystals, whereas, in our model the phonon transmittance is 
a strong function of frequency. 



VI. SUMMARY 

In summary, generic properties of elastic phonon transport at a disordered interface were 
studied by computing the exact scattering properties of a single disordered atomic layer 
sandwiched between identical waveguides. The results show that phonon transmittance is a 
strong function of frequency and the disorder correlation length. At frequencies lower than 
the van Hove singularity the transmittance at a given frequency increases as the correlation 
length decreases. At low frequencies, this is reflected by different power-laws for phonon 
conductance across correlated and uncorrelated disordered interfaces which are in approxi- 
mate agreement with perturbation theory of an elastic continuum ||. We also investigated 
simple mosaic and two-colour models of the interface and showed that phonon transmittance 
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through interfaces with correlated and uncorrelated disorder can be understood in terms of 
phonon conductances of structures with regular interfaces. 
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FIG. 1. A typical scattering geometry, where L(R) refer to sites on the faces of left(right) lead. 
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FIG. 2. Phonon transmittance for the (a)fcc model(^ = 5.65), (b)square scalar model(£ = 20), 
(c)estimates based on perturbation theory. In the (a) and (b) insets the corresponding normalised 
DOS versus uj is plotted. Frequency is measured in units of n/m. Inset (c) shows a single 
realization of rrij versus j for a line of correlated and uncorrelated masses(cr = 0.3 in all plots). 
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FIG. 3. Phonon transmittance for a = 0.3 and 0.1. Curves (1), (2) correspond to correlated 
(£ = 15) and uncorrected disorder for a = 0.3, whereas curves (3), (4) are for a = 0.1. 
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FIG. 4. Phonon transmittance for correlated disorder and regular models as explained in the 
text. In both (a) and (b) curves (1) and (3) correspond to correlated disorder for a = 0.3 and 0.1 
(£ = 15). In (a), curves (2) and (4) are for the regular interlayer with a' = 0.3 and 0.1. The results 
of the mosaic model are shown in (2b) and (4b). 



16 



1.0 r 

0.8 
0.6 
0.4 
0.2 




0.0 ' 1 1 ' — ' 

0.0 1.0 2.0 3.0 

CO 

FIG. 5. Phonon transmittance for uncorrelated disorder, curves (1) and (3) for a = 0.3 and 0.1, 
and the two-colour disorder model. 
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